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Abstract 

A method for calculating the relativistic path integral solution via sum over 
perturbation series is given. As an application the exact path integral solu- 
tion of the relativistic Aharonov-Bohm-Coulomb system is obtained by the 
method. Different from the earlier treatment based on the space-time trans- 
formation and infinite multiple-valued trasformation of Kustaanheimo-Stiefel 
in order to perform path integral, the method developed in this contribution 
involves only the explicit form of a simple Green's function and an explicit 

path integral is avoided. 
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Typeset using REVTeX 



-mail:d793314@phys.nthu.edu.tw 

1 



I. INTRODUCTION 



Based on the perturbation expansion of path integral formulation, Feynman firstly in- 
troduce his famous diagram technique to give a neat interpretation of the terms in the 
perturbation series and calculate the quantities of quantum electrodynamics order by order 
Over the last five decades, Feynman's method has been successfully applied to diverse 
areas of physics and achieved many accomplishments ||. Nevertheless, the exact result of 
summing the perturbation series is still the aim of seeking because of many physical effects 
in which non-perturbative exact result plays the pivot's role. In this contribution, a method 
for calculating the relativistic path integral is given in which the exact results only involve 
the computation of some kind of moments Q n over the Feynman measure and summing them 
in accordance with the Feynman-Kac type formula. So clear and neat is the method that 
it provides us not only with an alternative approach but a completely diverse viewpoint for 
treating physical problems. As an application, we apply the formula to calculate the path 
integral solution of the relativistic Aharonov-Bohm-Coulomb (A-B-C) system. It turns out 
that the method presented in this paper is neat due to the avoidance of space-time and 
(Kustaanheimo-Stiefel) K-S transformation in directly performming path integral |4[]. The 
A-B-C case can serves as the prototype for the treatment of arbitrary problems via summing 
the perturbation series. 



II. PATH INTEGRAL SOLUTION BY SUMMING THE PERTURBATION 

SERIES 

The starting point is the path integral representation for the Green's function of a rela- 
tivistic particle in external electromagnetic fields 0-§J: 

G(x b , Xa ;E) = ^-J™dS J Vp{\)<S>[p{\)\ J V D x(\) ex P {-A E {x,±] /h} p(0) (2.1) 

with the action 
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Ae [x, x] = / dX 

J Xa 



where S is defined as 



-x 2 (A) - z(e/c)A(x) • x(A) - p(A) — ^ + p (A) — - 



2p(A) v ; v / ; w v / rv y 2mc2 



(2.2) 



S= f Xb d\p(\), (2.3) 

A a 

in which p(A) is an arbitrary dimensionless fluctuating scale variable, p(0) is the terminal 
point of the function p(A), and $[p(A)] is some convenient gauge- fixing functional J§-||. The 
only condition on <£>[p(A)] is that 

^p(A)$[p(A)] = l. (2.4) 

h/mc is the well-known Compton wave length of a particle of mass m, A(x) and V(x) stand 
for the vector and scalar potential of the systems, respectively. E is the system energy, and 
x is the spatial part of the (D + 1) vector x M = (x, r). 

The functional integral for x in representation of Eq. ( [2.1| ) can be interpreted as the 
expectation value of the real functional exp j— i <iA/?p(A)V(x(A))| over the measure 

K (x 6 ,x a ;A 6 -A a ) = j V D x(\)e nJ ^ l 2pW clMj?U ^J, (2.5) 

and the entire Green's function reduces to the following formula 



jf, poo r 

G(X6, x a ; E) = - — / dS Vp(\)$ [p(A)] e 
2mc Jo J 



x (expl-l [ Xb d\Pp(\)V(x(\)))\ p(0) (2.6) 



H JXa 

in which £ = (m 2 c 4 — E 2 )/2mc 2 , (3 = E/mc 2 with the notation (★) standing for the expecta- 
tion value of the moment * over the measure i^ ( x fe ; x a! A& — A a ). Eq. ( |2.6| ) forms the basis for 
studying the relativistic potential problems by the Feynman-Kac type formula. Although 
we have chosen the term V^(x(A)) to expansion, it has the aesthetic appeal on choosing 
convenient one according to which the most suitable term is expanded for calculation. 



Expanding the potential V(x) in Eq. ( |2.6|) into a power series and interchanging the 
order of integration and summation, we have 

ih /•-.„/■ — i J^b d\p(X)£ 



G(x fe , x a ; E) = —j o dSj Vp<t> [p] e 



x t { -^P- n ( ( T d\p(\)V(x(\))) ^ p(0). 



(2.7) 



n=0 '*" \ V Xa / / 

We see that the calculation of path integral now turns into the computation of the ex- 
pectation value of moments Q n (Q = dXpV(x)) over the Feynman measure and 
summing them in accordance with the Feynman-Kac type formula. Ordering the A as 
Ai < A2 < ■ • • < X n < \b and denoting x(Aj) = Xj, the perturbation series in Eq. ( |2.7| ) 
explicitly turns into 



n=0 



A"o(x 6 ,x a ; A fe - A a ) 



+ V — r / dX n c/A n _i • • • / d\i 
~ x \ nj Jx a Jx a Jx a 



II ^o( x i+i> x jS A J+ i — Xj) 

j=0 



Y[piV(-Xi)dxi, 
(2.8) 



i=l 



where A = A a , A n+ i = X b , x n+i = x 6 , and x = x . 

As an application of Eq. ( [2.7|) , let's apply it to the relativistic A-B-C system in three 
dimensions. In this case, we have the vector and scalar potentials 



A(x)=24 yf + fy } V (r) = - e - 
x 2 + y 2 r 



(2.9) 



where e x>y stands for the unit vector along the x, y axis, respectively. The perturbative 
expansion in Eq. (|2T%|) becomes 



£ { -^~P- n ( i t dXp(X)z) > ^.(x,,. x„: A/, A„ ) 



n=0 



+ 



00 / q„2 " 



n=l 



> 2 

IT 



At r\ n 

dX n J dX n _i 

Aa J A a 



A 2 



Aa 



dX x 



II Kofej+i, Xj; Xj + i — Xj t 

j=0 



Up 

i=l 



<ix,- 



(2.10) 



The corresponding amplitude Kq takes the form 



f o -i f Xb d\ -™ x 2 (A)-i£A(x)-x(A)-p(A)#^si 

K (x b ,x a ,\ b -\ a ) = JV 3 xe nJ ^ N A) c (2 n) 

where a = e 2 /Tic is the fine structure constant. We now choose $ [p] = <5 [p — 1] to fix the 
value of p(A) to unity. The path integral in Eq. ( |2.7| ) becomes 

G(x fe , x a ; £) = — — / dSe~T> s { K (x b , x a ; S) 
2Mc Jo 



n=l 



II ^o(xj+i, Xj-; Aj+i - A,-) 

3=0 



TT *^Xj 

1 = 1 1 



(2.12) 



We observe that the integration over S is a Laplace transformation. Because of the convo- 
lution property of the Laplace transformation, we obtain 

ih 



G(x b ,x a ;E) 



2mc 



x 



Go(x6,x tt ;£) + £ (x) / 



Il G! o(x i +i ) x J ;f) 

3=0 



i r C?Xj 
1=1 1 



(2.13) 



with Go(x6, x a ; £) is the Laplace transformation of pseudopropagator K (yL b , x a ; A& — A a ). 
Let's first analyze the influence of A-B effect on the Go(x b , x a ; £). Introducing the azimuthal 
angle around the A-B tube 



<p(x) = arctan(y/x) 



(2.14) 



the components of the vector potential can be expressed as 



Ai = 2gdi<p(x) 



[2.15) 



The associated magnetic field lines are confined to an infinitely thin tube along the z-axis: 



B 3 = 2g€ 3ij d i d j (p(x) = 47t#5(x_l), 



(2.16) 



where xj_ stands for the transverse vector xj_ = (x,y). Note that the derivatives in front of 
<£>(x) commute everywhere, except at the origin where Stokes' theorem yields 

J d 2 x (d x d y — d y d x ) (p(x) = j> d(p = 2tt. (2.17) 

The magnetic flux through the tube is defined by the integral 

Q = f d 2 xB 3 . (2.18) 



This shows that the coupling constant g is related to the magnetic flux by 



Inserting = 2gdi(p(x) into the action of Eq. ( p. 11 ), the magnetic interaction takes the 
form 

Anag = ih{3 f S dXif(X), (2.20) 

J 

where y?(A) = </?(x(A)), (p = dip/dX, and (3$ is the dimensionless number 

A = - 2 -r-- (2.2i) 

The minus sign is a matter of convention. Since the particle orbits are present at all times, 
their worldlines in spacetime can be considered as being closed at infinity, and the integral 

l r s 

k = — dXip{X) (2.22) 
2n Jo 

is the topological invariant with integer values of the winding number k. The magnetic 
interaction is therefore purely topological, its value being 

A mag = ihf3 2kn. (2.23) 

The influence of A-B effect in the Green's function Go(x{>, x a ; S) is as follows. In the lacking 
of A-B effect, the Green's function 

G (x i+1 ,x,;£) = E gf\r^T J -£)Y lk ^ +1 )Y l l{± 3 ) 1 (2.24) 

n { r j+i r j) 1 l =0k =-l 



where Yzfc(x) is the 3-dimensional spherical harmonics Yzfc(x) and the is the radial Green's 
function of a particle moving in a centrifugal potential given by ]7],|8| 



p -m(r] +1+r ])/2ns T fm r j+1 rj \ 



(2.25) 



The notation / denotes the modified Bessel function. With the following formulas [p. 166, 

p.210, P 2i2, m, 



P i Q,/3) (*) = SrSf^ a + + n+l;a + l;±=* 
r(i+n) pf-uy^ _ r(i+n+p /s-lV p(M)/„\ 

r(l+n-0 " W — r(l+n) V 2 J ^n-J W 



(2.26) 



^ r(i+n-0 n v-^y — r(n 
where P£(x), Pjf*'^(x) are the associated Legendre polynomial and Jacobi function and F 
the hypergeometric function, it is not difficult to prove the following result 



(cos^/2sin^/2) K P/_X ; (cosi 



(2.27) 



The angular part of Eq. (|2.24 ) turns into 



E 



2/ + ir(i + z- fc) r (i + z + /c) 



47T 



r 2 (i + z) 



(cos Qj+i/2 cos9j/2 sm9j + i/2sm9j/2) 1 



(2.28) 



To go further, let's change the variable Z by defining I — k = q into g. It is easily to fine that 
the Green's function of Eq. (2.24) becomes 



/;/ 



2 (q + k) + 1 T (1 + q) T (1 + q + 2k) 



An 



T 2 (l + q + k) 



JKfj+i-'Pj) 



x (cos 6 j+1 /2 cos dj/2 sin j+1 /2 sin Q^lf P 9 (fc ' fc) (cos # i+1 )Pf ' fc) (cos d- t 



(2.29) 



7 



with gty k being the radial Green's function 



/o s v^+w-^U 5 J- 

Let us invoke the Poisson's summation formula [p. 469 ||] 



OO oo 

£ /(*:) = dy ^ e 2 ^f(y). 

J—oo „ 



k=— oo 



The entire Green's function Go(x&, x a ; £) containing the A-B effect becomes 



oo „ oo 
n f- 



fc=— OO 



(2.30) 



(2.31) 



x 



2(g + ^) + ir(l + g)r(l + g + 2z) 



47T 



T 2 (l + q + z) 



X 



(cos j+ i/2 cos 0^/2 sin i+1 /2 sin 0,-/2)* P (2 ' 2) (cos 6 j+1 )P^ z) (cos 0,). (2.32) 



The sum over all k in Eq. ( |2.32| ) forces z to be equal to fio modulo an arbitrary integral 
number leading to 



Co(x J+ i, Xj-; 5) 



^( r i+i r i) 1/2 



H H 5 , g+ i fc+/ 3 |( r j+i ) r i;^) 



q=0 k=—oo 



X 



2 (g + \k + Pq\) + 1 T (1 + q) T (1 + q + 2 \k + /3q\) 



An 



r 2 (l + g+|fc + /3 |) 



x (cos / 2 cos 0, / 2 sin J+1 / 2 sin 9j /2) 



|fc+/3o| 



X 



(2.33) 



with 



(0) / cn 



°° ^ P -f s -m(rJ +1 +»f)/2ftSj (m r j+1 rj 

5 V[ 2 (9+l fc +/3ol)+l] 2 -4a 2 /2 ^ 5 



(2.34) 



Using the orthogonality relations of Jacobi polynomials [p212, [|]], 

1 dx (1 - x) a (1 + xf P<*® (x) P£'V (x) 



2 Q+/3+1 T(a + n + l)T(/3 + n + 1) 



a + (3 + 2n + l n\T (a + (3 + n + 1) 
we perform the intermediate angular part of Eq. Q2.13D , it yields 

G(x 6 , x a ; = ^— 51 £ G 3i|fc+A) [(r 6 ,r a ;£) 



(2.35) 



q=0 fc=— oo 



X 



2 (g + |fc + A,|) + 1 T (1 + 9 ) r (1 + q + 2 I A; + (3 \) 



47T 



r2(i + g+|fc + /3 |) 



X 



^cos 6 b /2 cos # a /2 sin 6 b /2 sin # a /2) |fc+A)l p(l fe +A>U*+/%l) ( cos e ^ P (\k+M,\k +M ) (cos 0^ (2 .36) 



The pure radial amplitude G ni |fc + ^ | (rj,, r a ; £) has the form 

m 1 



Gq,|fc+/3 |( r 6> r a] 



ft {nr a ) 1/2 



^ ( mpe 2 \ n (n) 

2^ I ^2 ) # g+ |fc + /3 |l r &> r a,<>J 



n=0 



(2.37) 



with g^en by 



oo /*oo 



JO 



To obtain the explicit result of g^ k+f3o ^, we note that [|I 

dS 



n VHfc+/3o|( r i+ 1 > r i ; £) 

j=0 



(2.38) 



4 = 1 



o 5 



p \h S ) 



2 ^ efeJ_e-^ ooU, *J 2p f 2/tv f^ 
jo sinhz y sinhz 



(2.39) 



with /t = V m2 c 4 — E 2 /Tic. With the help of the integral formula |J, 
J o 2 



(2.40) 



we obtain the result 



^+\k +l 3 \(n, r a ; S) = I 9^ lk+M {n, r; S)g^ lk+M {r, r a ; £)dr 



(0) 



2^ /"OO 

— / zh(z)dz, 

K JO 



(2.41) 



where the function h(z) is defined as 



h(z) 



sinh 2 



-K(r b +r a ) coth z j 



V[2(9+|fc+A)l)+i] 2 -4a 2 I sinh , 



(2.42) 



The expression for #g™ | fc+(3o |(r&, r a ; 5) can be obtained by induction with respect to n, and is 



given by 



2" +1 1 f°° 

9g+ lk+M {n,r a ;S) = —r— / z n h(z)dz. 



nl k Jo 

Inserting the expression in Eq. ( p. 37] ), we obtain 

m 2 



(2.43) 



G q ,\k+(3 \(n,r a ;£) 



h {rbr a ) 1/2 



POO 

x / dze 
Jo 



I 2mf3e 2 \ 



_ e -K(r b +r a ) coth z j 



' 2Ky/r b r a 

Smhz~ V[ 2 (9+|fc+/3o|)+l] 2 -4a2 ^ s j n h z 

The integration can be done by the formula [e.g. ch. 9 0] 

t 



(2.44) 



I dy4 
Jo si 



,2vy 



sinh?/ 



exp 



(Ca + Cfe)cothy 



sinh?/ 



r((i + / z) /2-1/) 



(2.45) 



where and W^ tU are the Whittaker functions and the range of validity is given by 

Cb > (a > 0, 

Re[{l + /j,) /2 - u) > 0, 
Re(t) > 0, | argt |< 7r. 

We complete the integration and obtain 

1 



G q ,\k+f3 \( r b, T a \ £) 



mc 



(r b r a ) ^m 2 c* - E 2 
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r 1/2+ y/[2(q + \k + A)|) + 1] - 4a 2 /2 - Ea/Vm 2 ^ - E 2 
r (l + V[2(g+|fc + /? |) + l] 2 -4a 2 ) 



/ 2 



x M ^a/^ct . /?^IT^T^Z7li „ r-Vm¥-£V s . (2.46) 



EQ/Vm^i-i?^, v /[2( Q +|A:+/3o|)+l] 2 -4 a 2/2 ^ TlC 

The entire solution of path integral becomes 

ih mc 



G(x b , x a ; E) 



2mc 4irr b r a ^/m 2 c 4 - E 2 



oo oo 
><E E 



r 1/2 + J[2(q +\k + + If - 4ct 2 /2 - Pa/V^c 4 - E 2 



9=0fc=-oo I r(l + v /[2(g+|fc + /3 |) + l] 2 -4a 2 



XM 



X 



^r(l + + g + 2|fc + /? |) [2 (g + |fc + A)|) + 1] 

r 2 (l + g+|A; + /?o|) 

x e 4fc( ^" Va) (cos # b /2 cos # a /2 sin # 6 /2 sin o /2) |fc+A)l 



x P,(l fc + /3 ''l'l fc+/3o l)(cos^)P ( J |fc+/3o| ' |fc+/3o|) (cos^). (2.47) 

This result is given in Refs. ||, and p. 304 || in the first time where the same result must 
invoke the complicate space-time and multi-valued K-S transformations to perform the path 
integral. In present paper, this procedures is avoided and can be applied to arbitrary po- 
tential problems. 
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III. CONCLUDING REMARKS 



In the paper, a method for calculating the relativistic path integral involved essen- 
tially the computation of the expectation value of convenient moments Q n , such as Q = 

d\p(X)V(x) if we expands the term in action potential term, over the Feynman measure 
and summing them in accordance with the Feynman-Kac type formula is given. As an real- 
ization, the path integral solution of relativistic A-B-C system is given. Different from the 
former treatment in Ref. M, where the same problem must invoke the complicated space- 
time and the multi-valued K-S transformations to perform path integral, the merits of the 
method used in the paper is that it involves only the explicit form of some known Green's 
function and explicit path integral is avoided. The A-B-C system can serves as the prototype 
for the treatment of arbitrary problems via summing the perturbation series. It is our hope 
that our studies would help to achieve the ultimate goal of obtaining a comprehensive and 
complete solutions in perturbation series based on the path integral of quantum mechanics 
and quantum field theory, including quantum gravity and cosmology. 
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